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ABSTRACT 

In this study, we report on the evolution of particle size distributions (PSDs) and habits 
as measured during slow, Lagrangian-type spiral descents through deep subtropical and 
tropical cloud layers in Florida, Brazil, and Kwajalein, Marshall Islands, most of which 
were precipitating. The objective of the flight patterns was to learn more about how the 
PSDs evolved in the vertical and to obtain i n formation of the vertical structure of micro- 
physical properties. New instrumentation yielding better information on the concentra- 
tions of particles in the size (D) range between 0.2 and 2 cm, as well as improved particle 
imagery, produced more comprehensive observations for tropical stratiform precipitation 
regions and anvils than have been available previously. Collocated radar observations pro- 
vided additional information on the vertical structure of the cloud layers sampled. 

Most of the spirals began at cloud top, with temperatures (T) as low as — 50°C, and 
ended at cloud base or below the melting layer (ML). The PSDs broadened from cloud top 
towards cloud base, with the largest particles increasing in size from several millimeters 
at cloud top to one centimeter or larger towards cloud base. Some continued growth was 
noted in the upper part of the ML. Concentrations of particles less than 1 mm in size de- 
creased with decreasing height. The result was a consistent change in the PSDs in the ver- 
tical. Similarly, systematic changes in the size dependence of the particle cross-sectional 
area was noted with decreasing height. Aggregation — as ascertained from both the changes 
in the PSDs and evolution of particle habits as observed in high detail with the cloud par- 
ticle imager (CPI) probe — was responsible for these trends. 

The PSDs were generally well-represented by gamma distributions of the form N — NorD^e~^^^ 
that were fitted to the PSDs over 1 — km horizontal intervals throughout the spirals. The 
intercept (iVor ), slope (Ap ), and dispersion {fx) values were derived for each PSD. Expo- 
nential curves {N — = 0) were also fitted to the distributions. The Ap val- 

ues for given spirals varied systematically with temperature as did the values of A (expo- 
nential), and the data generally conformed to values found in previous studies involving 
exponential fits to size distributions in mid-latitude frontal and cirrus layers. Consider- 
able variability often noted in the PSD properties dining the loops of individual spirals 
was manifested primarily in large changes in iVop and Nq, but fx, Ap and A remained fairly 
stable. Temperature is not found to be the sole factor controlling Ap or A but is a primary 
one. Direct relationships were found between Ap and ATop or Ap and fx for the gamma dis- 
tributions and A and A^o for the exponential. The latter relationship was not found as dis- 
tinctly in earlier studies; observed PSDs in this study had better fidelity with less scatter. 

The p values changed monotonlcally with T over the range of temperatures and were di- 
rectly related to A^or or Ap , thereby reducing the number of variables in the PSD func- 
tional equation to two. In the upper part of the ML, Nq and A continued to decrease, and 
in the lower part these values began to increase as the largest particles melted. 

We developed general expressions relating various bulk microphysical, radar, and radia- 
tive transfer — related variables to A^or and Ap , useful for both tropical and mid-latitude 
clouds. These relationships facilitate the specification of a number of bulk properties in 
cloud and climate models. The results presented in this paper apply best to temperatures 
between 0 and — 40°C, for which the measured radar reflectivities fall in the range of 0 to 
25 dBZe. 
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1. INTRODUCTION 


During 1998 and 1999, four field campaigns were conducted by the Tropical Rainfall 
Measuring Mission (TRMM) to evaluate the performance of the TRMM radar and ra- 
diometer retrieval algorithms. These field programs also provided validation data for TRMM 
mesoscale and regional — scale models. The experiments were conducted in snbtropical 
and tropical regions: Texas and Florida (TExas-FLorida UNderfiights, TEFLUN-A and 
B), Brazil (Large Scale Biosphere- Atmosphere Experiment, LBA) and Kwajalein, Mar- 
shall Islands (Kwajalein Experiment, KWAJEX). Measurements were acquired using mnl- 
tipolarization gronnd-based Doppler radars, rain gauges, and in — situ and overflying air- 
craft. As part of the validation effort, in TEFLUN-B, LBA, and KWAJEX, the University 
of North Dakota (UND) Citation aircraft employed state — of — the-art instrumentation to 
acquire in — situ measurements. In this study, we report on microphysical data acquired 
by the Citation in deep tropical anvils and stratiform precipitating clouds and develop pa- 
rameterizations for use by the remote — sensing and modeling commnnities. 

Past in — situ microphysical observations in tropical ice clonds over the range — 80 < 

T < -20° C have provided indications of total particle concentrations {Nt) and ice water 
contents (IWC). In the primary studies, listed in Table 1, particles with minimum diam- 
eter of between 10 and 100 pm up to about 1 mm were measured. Griffith et al. (1980) 
found that the IWCs in three cirrus cloud decks sampled at multiple levels during the GARP 
Atlantic Tropical Experiment (GATE) generally increased from near cloud top (11.5 to 
13.0 km) to near cloud base (7 to 9.5 km); peak, level-averaged IWCs reached 0.05 g m“^. 
Knollenberg et al. (1982, 1993) reported IWCs as high as several hundredths g m“^ in anvil 
cirrus near Panama and in cirrus at the top of a cyclone off the northwest coast of Aus- 
tralia, both at r — 80°C. Takahashi and Kuhara (1993) used videosondes to character- 

ize the properties of cumulonimbus clouds over Pohnpei, Micronesia at T values as low as 
-80°C. They found that peak IWC and maximum diameter decreased with height (T) and 
approached the amonnts reported by Griffith et al. (1980) and Knollenberg et al. (1982) at 
lower values for T. Pueschel et al. (1995) reported on ice particle measurements in a ty- 
phoon over the range — 40 < T < 0°C. The highest concentrations and smallest particles 
were observed near -40°C. The microphysical characteristics of convectively generated ice 
clouds over the range — 70 to — 20°C, from observations near Kwajalein, and dnring the 
Central Eqnatorial Pacific Experiment (CEPEX), were described in Heymsfield and Mc- 
Farquhar (1996) and McFarquhar and Heymsfield (1996). Average IWCs over this range 
of T increased from about five thousandths to several tenths gram m“^, and IWC, cross- 
sectional area, and median mass-weighted diameter increased downward below cloud top. 
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Stith et al. (2002), from the TRMM field campaigns described below, showed that there 
was a wide variety of particle types observed both in the convective and stratiform regions 
of tropical ice clouds at temperatures from 0 to about -20° C, using data primarily from 
the cloud particle imager and 2D-C probes. Aggregation was observed to be a primary 
growth process, as was also the case in the Heymsfield and McFarquhar (f996) study. 

Although remote — sensing and modehng applications require knowledge of the PSDs in 
tropical clouds, relatively little data have been available to develop parameterizations for 
these PSDs. The PSDs of cirrus anvils and detached cirrus sampled during CEPEX were 
characterized by McFarquhar and Heymsfield (1997). They parameterized concentration 
in terms of the melted equivalent diameter {D^) as a function of IWC and T. A limitation 
of the data set was that concentrations of particles with size above approximately 0.1 cm 
were generally below the detection li m it of the 2D-C probe. 

This discussion of previous studies implies important gaps in our knowledge of tropi- 
cal ice cloud microphysical properties. In particular, few observations have been made of 
PSDs for particle sizes above 0.1 to 0.2 cm, or for T in the range 0 to — 20°C. There are 
also few observations of how PSDs vary in the vertical — other than those made by combin- 
ing horizontal flight legs at different altitudes — or of microphysical properties as related to 
radar echoes. Furthermore, little information is available on particle habits. The data in 
this study fill some of these gaps. Spiral descents were used to characterize microphysical 
properties in the vertical for temperatmres in the range — 50 to 0°C. Recently developed 
probes that provide high-quality particle habit information, especially for the smaller par- 
ticle sizes, were used. Also, considerably larger samphng volumes in the larger sizes pro- 
vided more reliable measurements of PSDs than those previously obtained for tropical ice 
clouds. 

In section 2, we describe new instrumentation used for this study. We present our ob- 
servations in section 3. IWCs and other bulk properties derived from the PSDs are given 
in section 4; we also develop parameterizations to represent these properties. In section 5, 
we compare our results to earlier observations in mid-latitude and tropical ice clouds. We 
conclude in section 6 by summarizing the results of this study. A list of symbols used for 
the different cloud microphysics variables is given in Appendix A, and the methods used 
to calculate particle mass and terminal velocity are given in Appendix B. 

2.INSTRUMENTATION AND DATA PROCESSING 

The microphysical data sets collected during the TRMM field campaigns probably con- 
stitute the most complete set of in — situ data in subtropical and tropical clouds to date. 
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On the UND Citation aircraft, particle shapes were imaged and size distribntions mea- 
snred from about 20 pm to >2 cm. 

A Stratton Park Engineering Co. (SPEC), Inc., cloud particle imager (CPI) probe pro- 
vided detailed information on the shapes and sizes of particles from approximately 20 pm 
to above 1 mm, with 2.3 pm resolution. Because the CPI sampling volume is relatively 
small and is still the subject of study, the size distributions obtained from the CPI are not 
used here. 

PSDs, along with low — resolution particle imagery, were obtained from the Particle Mea- 
suring Systems (PMS) 2D-C and the SPEC high volume particle spectrometer (HVPS) 
probes. The PSDs were measured from 33 pm to above 1 mm in 33 pm increments by the 
2D-C, and from 0.2 mm to about 5 cm in increments of 0.2 mm from the HVPS, with com- 
posite PSDs generated from the two probes by finding a point at which the size distribu- 
tions from each overlapped, or where the concentrations were comparable. The overlap 
size usually occurred from 1 to 2 mm. The nominal 2D-C probe sample volume is given 
by the probe’s array width (about 1 mm) times the separation between the probe’s arms 
(6.1 cm) times the Citation’s true airspeed (about 120 m s“^), or about 7 1 s“^. The “par- 
ticle reconstruction” technique was used to extend this sample volume by an amount that 
increased roughly linearly with size (see sample volume values in Heymsfield and Parrish, 
1979), by considering particles that were partially outside of the 2D-C probes’ sample vol- 
ume. The HVPS swept out a sample volume given by the array width (about 5 cm) times 
the separation between the probes’ arms of about 20 cm, times the aircraft’s true airspeed, 
or about 1 m“^ s“^. 

In general, the 2D-C probe worked well during the various TRMM field campaigns, al- 
though concentrations between 33 and 100 /xm are generally considered to be inaccurate 
because of questions related to the probe’s sample volume in this size range. The HVPS 
worked poorly in the TEFLUN experiments, worked intermittently in LBA, and worked 
well in KWAJEX. Therefore, examination of the vertical variability of the PSD empha- 
sizes the KWAJEX spirals. Useful CPI imagery was obtained from the spirals in all held 
programs. 

Key imaging probe data required for our analysis included the maximum particle diam- 
eter (D), cross-sectional area (A), area ratio (A^ = concentration per size bin 

4 

(Ni). The diameter was derived as the “true” maximum diameter — not the maximnm di- 
ameter along the array or Eight direction axes, as in previous calculations. 

State parameters were obtained from the aircraft’s standard suite of instruments. The 
T measurement, from a heated Rosemount sensor, is accurate to approximately ± 0.5°C. 
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A more complete discussion of the Citation state parameter instrnmentation appears in 
Stith et al. (2002). 

We used Doppler radar data from the ER-2 and DC-8 aircraft to characterize the reflec- 
tivity structure of the cloud layers during the Citation sampling, and ground-based S band 
(10.7 cm) radar data to examine temporal evolntion of the cloud layers during the spirals. 
The ER-2 Doppler radar (EDOP) transmits at a wavelength of 3.1 cm (G. Heymsfield, 

1996), whereas the airborne mapping Doppler radar (ARMAR) on the DC-8 transmits 
at a wavelength of 2.2 cm (Durden et ah, 1994). The vertical resolntion of these radars 
is 80 m (ARMAR) and 75 m (EDOP); the m in i m um detectable reflectivity as used here 
is approximately — 10 dBZe; and the precision of the velocity estimate (particle fall ve- 
locity + air velocity) is 0.1 m s“^. The vertical motions of each aircraft have been sub- 
tracted from the velocity estimates by using the Doppler velocity information in the low- 
est (surface) rangegate, together with the aircraft navigational system data to remove the 
aircraft vertical motion. Some random components of vertical motion conld not be sub- 
tracted but were eliminated through averaging ^ , which resulted in an absolnte accuracy 
of the Doppler velocities of ± 0.3 m s“^ for the ARMAR radar but was 0. 5-1.0 m s“^ for 
EDOP. The accuracy of the reflectivity measurements from each aircraft radar was ap- 
proximately ± 1 dB. No aircraft Doppler data were available for three of the flights ex- 
amined later but data from the ground-based radar for these flights was available. The 
ground-based radars included the NCAR S-pol radar during TEFLUN-B and LBA and 
a similar S-band radar during KWAJEX. (The characteristics of the KWAJEX radar are 
presented in Schumacher and Houze, 2001). 

3. OBSERVATIONS 

A. Sampling Strategy and Limitations 

The data reported here were collected primarily in anvils and dissipating thick strati- 
form precipitation regions. As pointed out by Field (1999), “it is quite common for an air- 
craft run at a given height to encounter a cloud that thins, breaks, and perhaps thickens 
again. If microphysical data are averaged over such a traverse it becomes difficult to disen- 
tangle the effects of differing local conditions on the evolution of the particle popnlation” . 

It is therefore difficult to assess how the properties of a cloud change in the vertical from 
data averaged for different levels. 

To examine vertical variability but particularly particle evolution in the vertical, each 
case stndied here involved a Lo and Passarelh (1982) type Lagrangian spiral descent through 

^ Possible upward motions in and above the vicinity of the radar bright band of order ten 
cm s“^ wonld not have been eliminated by averaging but will not significantly inflnence 
the resnlts. 
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a cloud to observe the evolution of the PSDs. A vertical profile using this technique starts 
aloft in a horizontally extensive ice clond area, and an aircraft is placed in a constant bank 
angle at a constant descent rate of about 1 m s“^. The aircraft spirals downward at ap- 
proximately the mean fallspeed of the snow and the loops of the spiral drift with the wind. 
Ideally, if conditions are quasi-steady and the properties of the atmosphere are fairly uni- 
form over a length scale somewhat larger than the diameter of the loops, then the aircraft 
largely samples particle evolution from the height change in the size distribution proper- 
ties. The analysis of particle size spectra can be performed by averaging spectra over a 
complete loop of the spiral. This serves to average any horizontal inhomogeneities. An- 
other approach is to compare particle size spectra at various heights that occurred in the 
same aircraft heading or sector of different loops. Ideally, if all the loops are the same size 
and the aircraft descends at about 1 m s~^, then each point of the aircraft trajectory cor- 
responds to the trajectory of a 1 m s~^ particle. 

Horizontal gradients in wind velocity, wind shear and dispersion of ice particle fallspeed 
creates a difficult sampling problem under the best of circumstances. For example, even if 
the cloud dynamics and microphysics are in a quasi-steady state, there are usually mesoscale 
areas of precipitation embedded within widespread cloud layers which move and change 
with time. Our observations are much more prone to sampling errors and misinterpreta- 
tion than earlier studies that used Lagrangian spiral descents to sample stratiform cloud 
layers. This is in part because temporal evolution was much more pronounced — the du- 
rations of the spirals were comparable to the life cycles of the convective clouds that gen- 
erated the cloud layers, such that the upper portions of the cloud we sampled could have 
evolved substantially by the time we sampled the lower portions of the clouds. Further- 
more, horizontal inhomogeneities are possibly more pronounced than in widespread cloud 
layers-small convective elements could have contributed concentrated but small-scale re- 
gions embedded within broader cloud layers. 

The vertical variations in the PSD properties that are presented in this section therefore 
may not represent the properties at any given time. What is being observed is the changes 
downwards in the part of the PSD falling at the descent velocity of the aircraft, whether 
or not there is substantial evolntion of the cloud properties above the sampling level. For 
each case studied, the extent of evolution of the radar echoes over the course of the de- 
scents will be characterized and the impact discussed in Section 5. Note also that horizon- 
tal variability in the properties of the PSDs dming the course of each spiral may be more 
extensive than in earlier studies, a factor that will also be discussed in Sections 3C and 5. 
And lastly, the spiral descents occupied less than five hours of aircraft flight time whereas 
there is almost 100 hours of data collected during horizontal traverses. This broader set 
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of data will be used in the future to study the relationship between radar reflectivity and 
properties of the PSD for tropical ice clouds in forthcoming articles planned by a nnmber 
of investigators. 

B. Overview of the Cloud Layer Properties 

The cloud top and base heights and temperatures measured by the aircraft during the 
spirals, the cloud optical depths estimated from the microphysical probes, and the changes 
in the radar reflectivities of the regions of the clouds sampled during the spirals are pre- 
sented in Table 2. Most spirals commenced at cloud top, with four descents initiating at 
temperatnres below — 35°C and five above — 20°C (Table 2). In the coldest case, T — — 
50°C. The spirals ended at clond base, or in rain when the clond base was below the ML, 
and spanned a cloud depth that was usually between 3.0 and 4.5 km. When possible, air- 
craft descent rates were adjusted to match the approximate mean particle fallspeeds 

both in the ice (1 to 2 m s“^) and rain (5 to 7 m s“^) regions. These attempts were suc- 
cessful in KWAJEX but not in LB A (Table 2, 

Each loop of the spiral spanned a diameter of 5 to 10 km, with the number of loops vary- 
ing from 3 to 23, exceeding 10 for all KWAJEX cases. The Citation drifted several tens of 
kilometers with the wind during the sampling periods (Fig. 1), ranging from 6 to 45 min- 
utes with most in the 25 to 45 minute range. 

An average optical depth for each spiral descent was derived by integrating downwards 
from the top to the base of the spiral, twice the cross-sectional area of the particle popula- 
tion as measured by the 2D-C and HVPS probes (Table 2). The cloud layers sampled were 
close to or exceeded the value of 23 necessary to be considered within the “Deep Convec- 
tion” cloud category according to the International Satellite Cloud Climatology Project 
(ISCCP, Rossow and Schiffer, 1999). Om clouds did not fall within the optically thinner 
categories of “Cirrus” or “Cirrostratus” . 

As shown in Table 2, “Radar” from the ground-based S-band radars, nine of the ten spi- 
ral descents were associated with precipitation (deduced from radar) at some time during 
the spirals. Most of the clouds would have had smtace reflectivities above the 18 dBZe or 
so necessary to be detected by the TRMM radar. (TEFLUN-B data courtesy Ed Brandes, 
LB A data courtesy Walt Petersen, KWAJEX data courtesy Sandra Yuter). 

Coincident vertical profiles of the measmed radar reflectivities (dBZe) and reflectivity- 
weighted Doppler particle fallspeeds {Vz, positive downward) were available from over- 
flying aircraft during seven of the spirals. The measured Vz values are the sum of the air 
and particle velocities. The ARMAR or EDOP provided information dnring or shortly be- 
fore these spirals (Fig. 1). Note that the radar data were nearly instantaneous, although 
the Citation spirals required 30 minutes or so and involved aircraft drift. While a one-to- 
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one comparison of the radar and in-situ measurements is therefore not possible , pertinent 
information on cloud structure, measured dBZg, Vz, and presence of significant updrafts 
were obtained, at least for some portion of the spirals. Periods were identified when the 
Citation and DC-8 or ER-2 tracks overlapped (bold lines. Fig. 1), and radar data during 
and near these periods were examined to assess the adjacent cloud structure. 

Averaged vertical profiles of measured dBZg and Vz for various periods of aircraft co- 
incidence are shown in Fig. 2, left and right panels. The dBZg in the ice regions (4.5 km 
above mean sea level) fell in the 5 — 25 dBZg range and generally increased downwards. 

The trend differed for the 990819 case (top panels), in which the dBZg values were nearly 
constant with height. The Vz in ice were 1 to 2 m s“^, also usually showing an increase 
downward. Vertical air motions may have been present in the data from individual radar 
scans, but averaging over the relatively long horizontal intervals probably canceled out up- 
and downdrafts, with the exception of possible local upward motions of order 10 cm s~^ 
at and above the melting layer. Vertical air motions also contributed to the ± lu bounds 
of 50 to 100 cm s“^ shown in the figure. Lower Vz generally coincided with lower dBZg. 
Bright bands — e.g., melting layers — were noted in the data at approximately 4.5 km in 
five of the cases. Below 4.2 km in rain, the Zg values were 20 to 30 dBZg, except on 990911, 
when lower values for dBZg were measured. The values for Vz in rain were about 6 m s“^. 

To place some bounds on the extent of temporal evolution during the spirals, airborne 
and ground-based radar measurements were used to loosely examine the change in the 
mean radar reflectivity at the Citation aircraft mid-spiral height and location at several 
times during each spiral. We use the FDOP and S-pol data for the 980905 spiral, the S- 
pol radar data for 990217 and 990219, the ARMAR data for 990819 and 990830, and the 
Kwajalein radar data for 990822, 990823 and 990911. The changes in the mean radar re- 
flectivity values, shown in increments of 5 dBZg in the last column of Table 2, were gener- 
ally 5 dBZ or below. As shown in Section 4, changes in dBZg of about 5 dBZg correspond 
to changes in the ice water content or precipitation rate by a factor of about two. This 
change is modest and therefore we conclude that the spiral descents provide not only an 
indication of the evolution of the PSD in the vertical but an indication of the approximate 
properties in the vertical over the time periods of the spirals. 

There was also considerable horizontal inhomogeneity over the course of the various loops 
of the spirals, and the variability (±l<r bounds) about the mean values of measured dBZg 
in Fig. 2 (left panels) was also often 5 dB. 

C. Particle Size Distribution Properties in the Vertical 

In KWA JEX, the rate of descent of the Citation coincided roughly with Doppler radar- 
measured particle fallspeeds. The particles dominating the measured reflectivities were 
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presumably at the larger end of the PSDs and their evolution in the vertical was there- 
fore captured by the Lagrangian spiral descents. However, evolution in the small end of 
the PSDs was not captured although aspects can be inferred from changes in the size dis- 
tributions. 

PSDs averaged over various loops for the spiral of 990819 are plotted in Fig. 3. The 2D- 
C and HVPS data overlapped reasonably well at about 0.2 cm, thus providing nearly con- 
tinuous PSD measurements for particles from tens of microns in size to centimeters. Most 
particles were below 0.2 cm. For particles above 0.04 to 0.06 cm, the concentrations (plot- 
ted on logarithmic axes) decreased linearly with size (linear axes) , indicating that at sizes 
above 0.04 to 0.06 cm, the PSDs were exponential. Below 0.04 to 0.06 cm, the PSDs were 
linear when plotted on logarithmic axes, indicating that particle concentrations in this 
part of the distribution were distributed by a power-law. The maximum particle size in- 
creased downward until the ML, then decreased below as a result of melting. 

PSDs plotted as in Fig. 3 do not captme the horizontal and vertical variability observed 
during the spirals. We developed a different type of plot to capture these trends. The left 
panels in Figs. 4 and 5 show height/time cross-sections of dBZg along the tracks of the 
DC-8 (KWAJEX) and ER-2 (LBA) as these aircraft flew near to and over the Citation. 
The boxes in the panels depict approximate locations of these overpasses if they had occurred 
over a period of about two minutes , for reference purposes only, but in fact they occurred over 
a thirty or so minute period (Table 2). The right panels show color-coded renditions of 
the PSDs during the spirals. Composite PSDs, as in Fig. 3, were color-coded to denote 
concentration per unit diameter, N{D), with higher N in red and lower N shown in blue. 
Colors represent N{D) (the color key is at the right), and the D scale is along the abscissa. 
The spirals used in Figs. 4 and 5 were selected because they contained both reliable N{D) 
information and coincident radar data. An exception is the 990217 spiral (bottom panel. 
Fig. 5), where some loss of PSD information occmred because of a probe malfunction. 

Inspection of the right panels in Figs. 4 and 5 shows the presence of cyclical fluctua- 
tions, indicating horizontal variability. This conclusion is supported by the radar imagery 
in the left panels, which indicates horizontal variability within the Citation domain sam- 
pled, sporadic bright band signatures, embedded convection, and distinct fallstreaks. It is, 
therefore, not necessarily meaningful to infer particle growth from loop averages for PSDs, 
and there are questions related to whether horizontal inhomogeneity is atypical of tropi- 
cal stratiform precipitating regions and anvils and if the observations are representative of 
tropical cloud layers. 

Beginning at cloud top, particles as large as 2 mm (9909II) to 6 mm (990819, 990830) 
were observed. The sizes of the largest particles in the PSDs increased downward, a fea- 
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ture shown by the decrease downwards in the width of the dark bine color in the right pan- 
els of Figs. 4 and 5. Fig. 6 shows images from the HVPS probe for each of the 14 loops of 
the Citation as it spiraled from cloud top to base on 990822. The particles shown for the 
first few loops were at most a few millimeters in size and reach up to 1 cm in diameter dur- 
ing the lower loops. 

The diameters of the largest particles increased downward at 0.1 to 0.3 cm km“^, with 
lesser increases for the colder cases and greater increases for the warmer cases. As the Ci- 
tation descended at a rate similar to the fall velocity of the largest particles, the broaden- 
ing of the PSDs can be attributed to growth of the larger particles and not to size sorting. 
Although large particles could have been injected into the lower levels, giving the illusion 
of growth, the measured dBZg profiles adjacent to the spiral locations do not suggest the 
presence of nearby sources of large particles. However, on 990822, the proximity to con- 
vection may suggest local sources of particles. 

The concentrations of small ice particles, which dominated the total concentration Nt, 
were greater near cloud top but generally decreased appreciably downward in the cloud 
(Fig. 7). Aggregates, examples of which appear in the lower panels of Fig. 8, were almost 
certainly responsible for the preponderance of this decrease, although a portion of the de- 
crease may have been the result of sublimation and size sorting. 

The PSDs changed markedly through the ML (Figs. 4 and 5: 990819, 990830, and 990911 
spirals). The tops of the ML — marked with “OC, ” in these figures for the 990819, 990830, 
and 990911 spirals — occurred at about 4.5 km. The size of the largest particles, Dmaxj 
initially increased from about 1 to 1.5 cm, while the N of particles under 0.5 cm decreased 
significantly (see Figs. 4, 5, and 7). Aggregation in the upper parts of the ML accounts 
for these observations. ^ The bases of the ML, as defined by the height at which the sizes 
of the largest particles no longer changed appreciably, occurred at about 3°C (3.95 km), 
where the largest raindrops were 0.3 — 0.4 cm. 

D. Particle Habit Variations in the Vertical 

The CPI provided information on how the particle habits varied in the vertical and the 
relationship of these to the measured reflectivity structures shown in Figs. 4 and 5. This 
instrument is best suited for characterizing the habits of the particles smaller than 2 mm; 
larger particles are not detected as its sample volume is too small to sample the larger par- 
ticles. On 990819 (as shown in Fig. 8), numerous vapor-grown crystals were observed in 

^ Some decrease in the total concentration could have been caused by the smaller particles 
decreasing in size below the 2D-C detection threshold, by melting resulting in increased 
fall velocities and decreased concentrations or by evaporation of the small drops by pref- 
erential deposition of condensate on the larger ice particles. Calculations indicate that 
these were not significant factors. 


9 



the intermediate (400 to 600 /rm) and large (>800 /xm) sizes — inclnding columns, capped 
columns, hexagonal plates, and branched crystals — whereas the habits of the small (<100 
/xm) particles were not identifiable. Aggregates, with some riming in the larger sizes, were 
also observed, especially just above and near the top of the ML. No supercooled liquid wa- 
ter was detected. The presence of the bright band for this spiral (top left panel. Fig. 6) in- 
dicated that updrafts were weak (because strong updrafts would disrupt the bright band) . 
Weak updrafts are conducive to growth primarily through diffusion rather than riming. 
Radar echo — top heights of 8 to 9 km (from T — 20 to — 25° C) led to the initial formation 
of columns, then of capped columns as particles fell through the planar-crystal growth 
regime where T ranged from — 12 to — 18°C. In contrast, particles observed on 990822 
(shown in Fig. 9) were rimed in the intermediate and large sizes with some aggregates ev- 
ident, and these appeared more spherical in the small sizes. Measured radar reflectivities 
(Fig. 4, middle left panel) and Vz (not shown), from which vertical velocities above a few 
m s“^ could be assessed, indicated the presence of deep updrafts, which led both to ex- 
tensive riming and to complex crystal shapes often associated with the freezing of cloud 
droplets at low temperatures. 

The relationship between particle habits, the measured radar reflectivities, and the prox- 
imity to convection was also observed in the other cases. The deep anvil generated from 
nearby deep and extensive convection on 990217 (Fig. 5, lower left panel, and EDOP data, 
not shown) produced complex, rimed crystals and aggregates of rimed crystals (upper pan- 
els, Fig. 10). No supercooled liquid water was measured on 990217, thus, riming must 
have been acquired in the convective regions, then particles were advected into the anvil. 
Conversely, for the poorly organized and weak updrafts associated with the 990823 spiral 
(lower panel, Fig. 4), habits consisted of pristine cirrus — type crystals (e.g., unrimed side 
planes and bullet rosettes), reflecting the low cloud — top T, and capped columns, reflect- 
ing additional growth around — 15°C. 

E. Derived moments of the size distributions 

Various bulk properties were computed from the PSDs, including the ice water content 
(IWC), precipitation rate (R), radar reflectivity factor (Z), radar reflectivity assuming 
equivalent water spheres dBZg, and mean mass and reflectivity- weighted particle termi- 
nal velocity {Vm, Vz), as in Heymsfield (1977). In this subsection, we present calculations 
of these parameters 

1 . Estimation of Ice Particle Density and Mass 

Calculation of each of these variables depends on a knowledge of the ice particle den- 
sity and mass as a function of D, and a knowledge of the terminal velocity (Vj), which are 
needed to calculate i?, Vz, and Vm- The Vt, in turn, depends on m, the ice particle cross- 
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sectional area (^4) normal to its fall direction, and the drag coefficient. Obtaining the area 
is straight-forward; A was measured directly by the 2D and HVPS probes, which were ori- 
ented to measure the particles’ horizontal cross-sections. 

Calculation of m is more problematic. Mass is obtained from the general relationship: 

m = 

where is the effective density (particle mass divided by the volume of a circumscribed 
sphere). No direct measurements of IWC were available from the TRMM observations to 
constrain the estimates of pe, and the ARMAR and EDOP measurements were not suffi- 
ciently collocated with those of the Citation to constrain the calculations from the radar 
reflectivity. Although Pe can be estimated for some particles with pristine habits, little di- 
rect information could be ascertained for the more complex, rimed, and aggregated parti- 
cles. A comprehensive study was used to estimate p^, encompassing a combination of cal- 
culations, ARMAR observations, and observations of ice particle terminal velocities at the 
surface; we describe this effort in Appendix B. The following relationship was developed 
from this analysis and from the work of Heymsfield et al. (2002), hereafter H02: 

Pe = MA,)”D«, (2) 

where A^. is the area ratio (,r/ 4 )D^ (number of shadowed 2D-C or HVPS pixels divided by 
the number of pixels in a circle with the same D), k — 0.04, n — 1.5, and a — -0.5 (cgs). 
The Ar provided some information about pe (see H02) and the D dependence accounted 
for observations that the pe for aggregates generally decreases with D. Although this re- 
lationship fit a number of observational data sets (described in Appendix B), the variabil- 
ity in Pe can nevertheless be large. We infer that using Eq. (2) with these coefficients pro- 
vides an accuracy of ±50% of the true IWC, as ascertained by comparing values for dBZe 
(oc m^) as calculated from the PSDs with those measured by the ARMAR on the same 
days during KWAJEX. An optimistic uncertainty in the mean pe of ± 25% leads to the 
following approximate uncertainties: IWC, ± 25%; R, ± 45%; dBZg, ± 10% (for positive 
dBZe); ^Zi i 12%; and Vmi i 15%. We have not accounted for Mie scattering effects in 
these estimates; this effect minimally impacts the calculations at the ARMAR and EDOP 
wavelengths (Robert Meneghini, private communication). 

2. Derived IWCs and other moments 

Vertical profiles of IWC, Z, and Vz as calculated from the PSDs for three spirals 
(shown in Fig. 11) were selected as representative of the findings for the ice regions of all 
spirals. The equations used to calculate these properties are presented and parameterized 


11 



in section 4C and listed in Table 3. The calculations are derived for temperatures 0°C and 
below to avoid ambiguities associated with particle mass in the melting layer. Our main 
findings were as follows: 

1 The magnitudes of the IWCs are mostly from 0.1 to 0.5 g m“^, and change by only about 
a factor of two or three with height, much less than is observed for synoptically-generated 
ice clouds over comparable vertical depths (see Section 5 and Fig. 18). The horizontal 
variability noted in the IWC values of less than 50% during the course of a spiral is com- 
parable to that observed in earlier Lagrangian spiral descents reported on by Lo and 
Passarelli (1982) and Field (1999). The magnitudes of the IWCs, and associated cal- 
culated values for dBZe, were comparable to those for mid-latitude anvils (Heymsfield 
1986). 

2 The precipitation rates changed by less than a factor of two over the course of the var- 
ious loops of each spiral and by less than a factor of three over the cloud depths, with 
the exception of the 990822 case. This result is good given that our goal was to sample 
the same population of particles as they fell from upper to lower cloud levels, implying a 
nearly constant R with height. 

3 The trends of calculated dBZg with height differed markedly from the radar observa- 
tions, shown in Fig. 2. This is not surprising given that the radar data were acquired 
over only part of the Citation track and for a very limited portion of each spiral. The 
calculated peak values of dBZe, however, were comparable to the radar measurements. 

Mie scattering effects should be included in cases where calculations and measurements 
of dBZe can be made directly. 

4 The Vjn were of the order 1 m s“^, very reasonable for the calculated IWCs and Rs. 

4. PROPERTIES AND PARAMETERIZATIONS OF THE PSDs 

The focus of this section is to develop general relationships for the PSDs, and between 
the PSDs and various microphysical and radar-measured properties for use in cloud-resolving 
models, climate models, and remote-sensing applications. 

A. Form of the Ice PSDs 

In this subsection, we use curves fitted to the PSDs to quantify how the PSDs from the 
Lagrangian spirals varied with temperature and height, and to quantitatively assess how 
horizontal variability influenced the PSD. Not all spirals were used in this analysis, be- 
cause a complete set of 2D and HVPS data was not always available. 

Based on the work of Kosarev and Mazin (1991), Mitchell (1991) and others, gamma 
distributions of the form 
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N = NorD^'e 


( 3 ) 


-XtD 


were fitted to the PSDs with diameters 66 fim or two 2D-C probe size bins and above for 
each 1 km of horizontal distance, or about 20 m in the vertical and 8-9 sec of flight. The 
intercept (iVor ), slope (Ar ), and dispersion (/r) values were derived for each PSD by match- 
ing three moments (Kozu and Nakamura 1991; Zhang et al. 2001). We have chosen to use 
the first, second, and sixth moments, as this set provided the best fit over the entire mea- 
sured particle size range. (Symbols in the discussion that follows are dehned in Appendix 
A.) 

The fit coefficients were derived as follows. The pth moment M{p) of the observed PSD 
is given by 


M(p) = / p(D)DPdD. 

JQ 


( 4 ) 


The parameter fi can be derived by finding the only real root of the quartic polynomial 


(1 - F)/ + (8 - 18F)iJ,^ + (24 - 119F)fj,'^ + (32 - 342F);U + (16 - 360F) = 0, (5) 

where F — other fit parameters are derived from 


and 


M{l)ip + 2) 
M(2) ’ 


( 6 ) 


A^or 


M(l)Ar^'‘+^^ 
P(^ + 2) 


( 7 ) 


This method produces a ht that is also in good agreement with other moments that are 
not used in the fitting routine, such as the third moment, used for computing IWC. Also 
htted were exponential curves (/r = 0) in Eq. (3)), for reference to earlier studies. Here- 
after, the ht coefficients for the exponential will be noted by A and Nq, and for the gam- 
mas by Ar, iVor, and p. 

The values of A and Ar obtained for the various spirals showed a correlation with T (Figs. 
12A and 12B) and with each other (Fig. 13A), although there was considerable variability. 
The A and Ar values were about the same (Fig. 13A), except below A and Ar less than 15 
cm“^ and above A and Ar greater than 70 cm“^. The solid line in the hgure (and the ht- 
ted equation) shows where the values deviate. 
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For T < — 15°C, the A and Ap values were usually greater than 20 cm with the excep- 
tion of the 990822 case, which differed, presumably because of the proximity of the sam- 
pling to convection and possibly because it was an anvil. For values of T between -5 and 
-f5°C, the A values decreased to between 12 and 20 cm“^ while the Ap were somewhat 
lower. As T increased toward 0°C, the A values were asymptotic at 8 to 10 cm“^, similar 
to findings in earlier studies (Houze et al. 1979; Lo and Passarelli 1982), which indicated 
that A values did not decrease below 8 to 10 cm“^. The Ap values below approximately 4 
cm“i resulted from poor curve fits, as ascertained from the correlation coefficients for the 
curve fits. 

For individual spirals, the A versus T values and the Ap versus T curves fell along nearly 
straight lines, with slopes or that were within a factor of 2 of each other over all 
cases. Within loops of individual spirals, the A or Ap values usually fluctuated by at most 
± 50%, far less than the changes noted over the comse of a spiral (see insert, Fig. 12B). 

The A values conformed quite closely to the A equation from Ryan (2000) for mid-latitude 
ice clouds that has been converted to an exponential form here, 

A (cm“^) = Co exp^*^*^f*^^\ (8) 

where the cq and ci values adapted from the coefficients in Ryan are are 12.1 cm"’^ and 
— 0.0564, respectively, and the curve is plotted in Fig. 12A. Coefficients co and ci were de- 
rived from fits to the median TRMM values of the distribution of A versus T and Ap ver- 
sus T shown in Figs. 12A and 12B, along with curves fitted to the ± \a values. (The data 
from the outlying 990822 case was omitted from the fitted data set, and two curves were 
fitted to the Ap versus T data). Cmves representing these equations are plotted in Figs. 

12A and 12B and are listed at the bottom of Table 3. In and below the ML, A and Ap val- 
ues again increased, in much the same way as was shown for A values by Houze et al. (1979). 

We can assess how accurately T would have predicted Ap by using the TRMM temper- 
atures in the spirals and comparing the predicted to measured Ap (Fig. 13B). As is shown 
in the figure, this approach to predicting Ap leads to considerable scatter, especially for 
the outlying case on 990822, where virtually all of the points deviate widely from the line 
in Fig. 13B — thus, T is not an ideal predictor of Ap (or A) values. A similar interpretation 
results if the Ryan (2000) relationship is used to estimate A. 

The Nq and A^or values also decreased systematically with increasing T (Figs. 12C and 
12D). However, there is no strong dependence of Nq or Nqy on T, which is why Ryan (2000) 
also found no clear relationship between Nq and T, only noting that “the parameter A^o 
is a less systematic function of T with local geographic variations being evident.” . Fur- 
thermore, there is considerably more scatter than was found for A and Ap especially over 
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individual loops, where Nq or Nor often varied by up to one or two magnitudes (insert, 

Fig. 12D). Since IWC is proportional to Nq and JVor (Section 4C), the variability of the 
Nq and iVgr represents IWC variability, and indicates that there are large fluctuations in 
IWC across the loops of the spirals. 

The correlation coefficients (r) for the exponential fits (Fig. 12E) were generally high for 
T values between -35 and -5°C but not outside of this range. For A values above 80 cm“^ 
or below 15 cm“^, the quality of the exponential cin*ve fits decreased. Nevertheless, the r 
value averaged over all exponential fits was 0.96, signifying that the exponential fits were 
quite good. The gamma distribution fits produced even better correlation coefficients (not 
shown) with an average r value of 0.98, and eliminated the problem areas at high and low 
A. 

At the larger values of Ar , the tended to have positive values (subexponential distri- 
butions), whereas at the lower values of Ar /r had negative values (Fig. 13C) (superexpo- 
nential distributions). Fluctuations in the fj, values during the loops of individual spirals 
were of order of 1.0 (insert. Fig. 12F), signifying that there was considerable variability in 
the shape of the PSD in small sizes. The fj, values were close to -1 in the range 25 < A < 

70 (Fig. 13C), which is why the correlation coefficients for the exponential were relatively 
high in this range. The exception was noted on 990822, where iJ, values fell below -1. 

The values of fj, were highly correlated with both Ar and A^or (Figs. 13C-D), especially 
below 0°C. The curve fits shown in these figmres and listed at the bottom of Table 3 can be 
used to eliminate the A^or term in Eq. (3) and reduce the number of unknown variables in 
Eq. (3) to two. The fits were derived with and without the data for the 990822 case, with 
little difference found. 

An indirect relationship was found between the maximum measured particle size, Dmax^ 
and A or Ap (Fig. 14). ^ At least two particles were required in a given HVPS size bin to 
find a Djnax- The D^ax at high A or Ar (low T) values were several millimeters and at low 
values were 1 cm or above. Curves fitted to the mean and standard deviation values are 
indicated in the figure and listed at the bottom of Table 3, with little difference found with 
and without the inclusion of the 990822 data. Extrapolation of the D^^ax curve in Fig. 

14A to 2 cm or above returns A values of only 6 or 7 cm“^. We speculate that the absence 
of A values of less than 8 or 9 cm“^ In earlier studies is due to a rarity of aggregates larger 
than 2 cm. 

^Decreasing A or Ap values and hence increasing spectral breadth will generally lead to 
an increase in Dmax- The measured Dj^ax also depends on the probe sampling volnme 
and No or A^or (which displace the PSDs upward or downward), and this relationship 
embodies these factors implicitly. 
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The Nq and A or iVgr and Ap values were highly correlated for individual spirals (Figs. 

15A and 15B) and tended to decrease from large values at the top of the spirals to small 
values in the upper parts of the ML (indicated by o in each panel). Broadening of the PSDs 
with distance below cloud top by aggregation and with a corresponding decrease in Nq or 
iVor by depletion of the smaller particles can account for this trend, although there are is- 
sues related to whether evolution of small particles was actually being observed during the 
Lagrangian spirals. The increase in D by diffusional growth is very slow for large parti- 
cles even if the environment was supersatmated with respect to ice and cannot account for 
these observations. 

The 0°C level was denoted by a marked change in the trend of the Nq versus A or iVor versus 
Ar points with T (Figs. 15A and 15B). Between 0 and 1.5°C, corresponding to the first 
several hundred meters, the A or Ar values decreased slightly, reflecting slight increases 
of Dmax and decreases in the total crystal concentration (Figs. 4 and 5; 990819, 990830, 
and 990911 cases). Broadening in the large end of the PSDs is consistent with previous 
observations by Yokoyama (1985) and Willis and Heymsfield (1989), who reported that 
aggregational growth continued several hundred meters into the ML. The decrease in Nq 
or A^or is only partially accounted for by ice particles melting to sizes below the 2D-C de- 
tection threshold because particles of sizes of 50 /um and below have high densities and do 
not change their D values much by melting. Most of the decrease is attributable to deple- 
tion by aggregation. 

The lower part of the ML was characterized by an increase in the Nq and A^or values 
(from the o to A symbols in Figs. 15A and 15B). These increases suggest that there is lit- 
tle breakup of partially melted individual aggregates to multiple particles in this region, as 
the total concentration, from Nt — (for an exponential, Sekhon and Srivastava 1970), 
remained essentially constant. The increases in the A or Ar values in the lower part of the 
ML resulted from a dominance of melting by the larger, lower — density aggregates. At 
this height in the ML, most of the smaller ice particles have melted, thus, their size is basi- 
cally unchanged. The diameters of the larger, relatively low — density aggregates changed 
signihcantly. Thus, the A and Nq or Ar and Afor values increased slightly. 

B. Parameterizations of the Ice Particle Shapes 

The area ratio is an important property in our estimates of ice particle density given 
by Eq. (2) and in parameterizing area-related properties (terminal velocity, extinction) 
of the PSDs. The Aj. is expressed as a power-law function of D\ Aj. — aD^. The coeffi- 
cients a and b can be derived through curve fits to the A^ versus D data from the imaging 
probes, collected over 1-km intervals dming the spirals. As shown in Figs. 16A to 16D, 
the a and b coefficients are a strong function of altitnde and temperature, and not sur- 
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prisingly tend towards spheres (a —1, fe = 0) in the melting layer. The height trends con- 
form to the tendencies fonnd for midlatitude cirrus clouds by Heymsfield and Miloshevich 
(2002), who attributed the changes in height or temperature to aggregation and sublima- 
tion. Points for the 990822 case, shown as grey symbols in Fig. 16, differ somewhat from 
those for the other cases. Given that the A values decrease downwards, it can be inferred 
that the height and temperature tendencies are consistent with changes in A. This resnlt 
is co nfir med in Fig. 16E, which shows the changes in the coefficient a with A. A nearly 
monotonic relationship is fonnd between b and A, or between a and b (Fig. 16F), with rela- 
tively little scatter. Cnrve fits between a and A, and a and b, are shown in Table 3, both to 
the mean values and standard deviations to indicate the extent of the scatter in the valnes. 

C. Parameterizations of the Moments of the PSDs 

The set of eqnations used to calculate the bulk properties of the particle population, in- 
clnding the IWC, R, Z, Vm, Vz, the extinction coefficient, e, and the effective radius, r^, 
are shown in Table (3), with the symbols identified in Appendix A. They nse as input the 
values from the relationship between the various properties of the gamma distributions 
and temperature reported in Section 4A; these equations are listed at the bottom of Ta- 
ble 3. The equations make use of PSDs represented by gamma distributions (Eq. 1) which 
provides an accurate representation of the TRMM PSDs, although they can readily use 
exponential if the value of fj, is taken to be 0. The equations listed under “Full Equations” 
in the table are general and can use many of the earlier measurements of the properties of 
exponential size distributions for ice particles reported later in Section 5. Particles below 
the 2D-C probe detection threshold are included implicitly through the use of the gamma 
distributions, but the accuracy in these sizes cannot be estimated. The resulting errors 
might have a major affect on the lower moments (primarily on total concentration and 
cross-sectional area) of the distributions but not on the higher moments (IWC or dBZe). 

The IWC can be expressed as 

IWC = -kNo (9) 

6 Jo 

using the coefficients a, k, and n identified in section 4 A and at the bottom of Table 3. 

The resulting solution to the integral of Eq. (9) is shown in column 2 of Table 3 for inte- 
gration from 0 to oo (complete gamma function) . We also integrated this equation over 
the range 0 to (incomplete gamma function) , with a generic solution shown at the 

bottom of the table. For given A^or and Ap values, IWCs returned from the complete and 
incomplete gamma functions differed insignificantly. A simplified relationship derived us- 
ing the values of a, k, and n are presented in coliunn 3. 
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The relative accuracy of the estimates of IWC using the exponential and gamma PSDs 
is shown in Figs. 17A and 17B. The numerators given in each of these panels were derived 
by taking the IWCs, as calculated from Eq. (9), using the fitted Nq and A or iVgr and Ap val- 
ues and the TRMM average Ay versus D relationship. The denominators are the IWCs 
derived directly from the PSDs and measnred Ay values. Using exponential size distribu- 
tions to derive IWC values leads to significant departures in the ratios from unity. Dra- 
matic improvements were noted throngh the nse of gamma distributions, yielding values 
generally close to unity. Therefore, the gamma distributions provide a much better repre- 
sentation of the PSDs than do the exponential. 

If IWC is known from measurements or prognosed in a general circulation model (GCM), 
the corresponding Nq and A or iVop and Ap valnes can be readily found. The relationship 
between A and T or Ap and T at the bottom of Table 3 can be used to prescribe A or Ap ; 
given an IWC, Nq or Nq^ values can then be obtained from the IWC relations shown in 
Table 3. This technique is illustrated in Fig. 17C: Ap values are derived from T (bottom. 
Table 3); the iVop values and p values are taken directly from the gamma fits; and then the 
IWC is calculated. The ratios of the calculated IWC values to the values taken directly 
from the PSDs may then be derived. In essence, this method tests for the error that can 
be expected in A^or when IWC and T are known. The grouping of points near the bot- 
tom of the panel actually represents data points and shows that for the 990822 case this 
method greatly underestimated the IWCs because the Ap values were overestimated. For 
points for the other spirals the scatter is large but the method is generally accurate to ±50%. 
Use of the exponential (not shown) produces larger errors because the population of small 
particles is inaccurately prescribed by the exponential. More accurate estimates of Nq or 
iVop are not possible without better understanding of the factors responsible for the vari- 
ability in A or Ap . 

The median mass diameter, Dynmi a property related to the distribution of ice mass with 
size, can be found by modifying the results presented in Mitchell (1991) and is a function 
of Ap, ji, and the coefficients in the exponent of D in Eq. (9). The equation is presented in 
Table 3. 

The radar reflectivity {Z) can be calculated from 

rD 

Z = / N{D)DjdD, (10) 

where Dm is the melted equivalent diameter, found using the expression for as given 
in Eq. (2). Equation (10), for Z, is analogous to Eq. (9) for IWC, with the results shown 
in Table 3. Wavelength-dependent, Mie scattering effects were not considered in the in- 
tegral evaluation. As with IWCs, the gannna distributions provide a better match to the 
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values calculated from the PSDs than do the exponential (see Figs. 17D and 17E). Given 
values for Z and T, iVo or A^or can be found in a similar manner to that described for cal- 
culating IWC. However, this method may lead to very large errors (Fig. 17F). The equiv- 
alent radar reflectivity dBZg listed in Table 3 is derived assuming equivalent water spheres 
(Smith, 1984). 

The precipitation rate (R) was calculated by incorporating the following H^-dependent 
Vt relationship from H02 inside the integral in Eq. (9): 

^ - 2G)^[(3+6(n-l)+a)?,,-l] ^ 

^Pa 

where v is the kinematic viscosity of air, pa is the density of air, and a f and b f are the co- 
efficients nsed in fall velocity relationships of the form introduced by Mitchell (1996) and 
recently modified by H02 to account for aggregates. To solve Eq. (11), only one set of 
and bf coefficients can be nsed at this time to represent all sizes. We carefully examined 
the use of appropriate a/ and bf coefficients, and selected those shown in Table 3. The in- 
tegrated equations for R appear in Table 3. 

Equations for Vz and Vm also follow the development used for Eqs. (9) and (11). The Vt 
given by Eq. (11) is included in the integral in Eq. (9) and normalized by the IWC to ob- 
tain Vm', Vt in Eq. (11) is incorporated in the integral in Eq. (10) and normalized by Z to 
obtain Vz- Solutions of these equations appears in Table 3. Note that these fall- velocity 
expressions do not depend directly on Nq but do depend on p and thus indirectly on Nq. 

Also note that different sets of the a/, 6/ coefficients are used for Vm and Vz resulting from 
different parts of the PSD contributing to each. The Vm equation predicts values that are 
accurate to ±20% (mostly slight overestimates) for Dmm above 0.03 cm (>90% of the points) 
and the Vm equation predicts values that are accmate to ± 20% for Dmm from 0.02 to 0.2 
cm (83% of the points). 

Given Vz and Z values, the Vz expression can be used to derive Ar , iVor , and p val- 
ues. With the use of Z or IWC, the full properties of the exponential PSDs are known. 
Vertically-pointing Doppler radar data, with some time averaging to remove small-scale 
vertical motions , can therefore be used to obtain direct measurements of Ap and Nq with- 
out the need to make assumptions about the temperature dependence of Ap . 

The total cross-sectional area of a population of particles (Ac) is given by 

^ N{D)ArD^dD, (12) 

^It may not be possible to remove vertical velocities of order 10 cm s“^ associated with 
the melting layer and synoptic systems by averaging over long periods of time 
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and the extinction coefficient in visible wavelengths (when particles are large compared 
to the wavelength) is calculated from e = 2Ac- Solutions for these integrals are listed in 
Table 3, where the a and b coefficients can be taken directly from either the relationship 
to A shown in the bottom of Table 3 or from the average values of each coefficient listed in 
the table. 

The effective radius (re), a parameter that characterizes the radiative size of a popula- 
tion of particles, can be written, according to Fu (1996), as 


y/3IWC 
3ft Ac 


(13) 


where pi — 0.91 g m“^. Solution of this equation derives from Eqs. (9) and (12) and is 
given in Table 3. For the TRMM spirals, the absence of measurements from small parti- 
cles is not thought to significantly affect the values. Heymsfield and McFarquhar (1996), 
who did have measurements of size spectra down to about 10pm in size, showed that the 
re at the warmer temperatures in tropical anvils were dominated by the larger particles. 


5.DISCUSSION 


Lagrangian spiral-type descents and ascents have been used previously to examine the 
evolution of PSDs in mid-latitude layered and frontal clouds by Passarelli (1978), Lo and 
Passarelli (1982), Gordon and Marwitz (1986), and Field (1999) (hereafter P78, LP82, 
GM86, and F99, respectively). In these studies, values for T at cloud top (-30°C) and base 
(2°C) were comparable to those of the TRMM observations. PSDs were measured using 
ID-P (P78, LP82); or 2D-C and 2D-P (GM86); or the 2D-C and HVPS probes (F99). Dif- 
fusional and aggregational growth predominated in each of the spirals. 

The parameters Nq and A of the exponential fitted to their PSDs by each of the authors 
are plotted in Fig. 18. More scatter was evident in the Nq versus A values than we found 
for the exponential in our spirals because the probe sample volumes and upper size detec- 
tion limits in all but the F99 study were significantly smaller than those used in our study. 
For reference, the T values at cloud top and base are shown in each panel, and the approx- 
imate envelope of the Nq versus A values from the TRMM observations are given in the 
lower right panel. Also plotted in two of the panels are the A values derived from the mea- 
sured T from the TRMM observations (bottom of Table 3). 

The following similarities and differences may be noted between these and the TRMM 
observations as shown in Fig. 15A: 

1. The A values generally fell in the range 10 to 50 cm“^, similar to the TRMM observa- 
tions over the same T range. 
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2. The trends in the Nq versus A values were si mi lar to those we observed. One of the GM86 
and two of the LP82 cases first showed an increase in the Nq value with a decrease in the 
value of A at cloud top. This pattern was attributed by the authors to particles grow- 
ing by diffusion into the size-detection threshold of the ID-P and 2D-C probes and to ice 
nucleation. 

3. The TRMM parameterization does a reasonably good job of predicting the value of A 
when A decreases with T, at T below 0°C, and when sublimating bases are not present. 

4. The Nq values for most but not all cases were comparable to the TRMM Nq. Above 
0°C, A^o values for the two GM86 cases and those from the TRMM observations corre- 
sponded closely. R is in this region that the probe sampling volume issues are probably 
minimized. Values for Nq for several of the LP82 and P78 cases were below those mea- 
sured during TRMM. We believe that this results from an underestimate of the concen- 
trations of crystals in the smaller size range of the ID-P, a conjecture confirmed by com- 
parisons in GM86 of Nq values from the 2D-C and 2D-P probes. As reported in GM86, 
the ratio A^o(2D — P)/A^o(2D — C) = 0.21, whereas the ratio A(2D — P)/(2D — C) = 0.72; 
thus, Nq changes appreciably but A doesn’t. 

5. The range of IWCs sampled from these various spirals was comparable to those from the 
TRMM spirals. The lines of constant IWC in Fig. 18 were derived from 

IWC(g m“^) = ^ Noyr(4+^))Xio ^ comes from Sekhon and Srivastava (1970), with 
y=0.0056 and r]— — 1.1 taken from Brown and Francis (1994). For a given constant value 
of IWC and a value of A, A^o values were derived and the lines of constant IWC were plot- 
ted. 

6. Most importantly, the gamma functions would have provided better estimates of the 
PSDs for A below 20 to 30 cm“^. 

Sekhon and Srivastava (1970) derived the following relationships for snowfall collected 
at the ground: 


Nq = 2.50xl0"^R“°-^'^ (cm"^), A = 22.9i?“°-^^ (cm“^), (14) 

where A^o and A are for the resulting melted PSDs. Although no direct relationship be- 
tween Nq and A was derived, this relationship can be estimated by eliminating R in Eq. 
(11), yielding 


Nq = 3.61xlO“^A^°^ (15) 

We could have developed a more accurate relationship than that given by Eq. (15) if the 
original Nq versus A data points were available. In Eig. 19, Eq. (15) is plotted over the 
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range of R reported in that study. Also shown in the fignre are the Nq versus A points for 
the 990822 spiral and the corresponding “melted” Nq versus A points. These were obtained 
by htting exponential curves to the concentration versus melted equivalent diameters (from 
Pe in Eq. [2]) for each PSD. 

Our melted Nq versus A points differed somewhat from points along the curve given by 
Eq. (15), possibly for the following reasons. At T above approximately — 5°C, the A val- 
ues for the unmelted and partially melted spectra remain relatively constant at 9 cm“^. 

This regime could not be discerned from the small number of points used in the Sekhon 
and Srivastava analyses. Furthermore, partially melted snow and/or raindrops were in- 
cluded in about half of the samples. Inspection of Fig. 11 shows that for a narrow range of 
T values between 0 and 3°C (for each spiral within this range, the R values were approxi- 
mately constant), the A values varied by a factor of 2-2.5. Thus, for partially melted snow, 
there is no unique R versus A relationship, and the data points within the region of melt- 
ing are suspect. 

McFarquhar and Heymsfield (1997) (hereafter MH97) recently developed a parameter- 
ization for the PSDs of tropical ice clouds (produced by outflows of deep convection) as 
functions of T and IWC, with stated ranges of validity -70 < T < -20 °C and for 10“^ < IWC < 
1 g m“^. The intention of McFarquhar and Heymsfield was to produce reliable measures 
of both IWC and cross-sectional area. The authors emphasized in this paper the small 
particles that often dominate cloud extinction and optical depth. The PSDs, which were 
measured from about 10 pm to a size larger than 1 mm, depending on the particle con- 
centration (Table 1), were converted to N{Dm) for Dm < 1000 pm through use of mass- 
diameter relationships. They were then represented as the sum of a first-order gamma 
function, describing ice crystals with Dm < 100 pm, and a lognormal function, describ- 
ing larger crystals. The percentage of the IWC comprising small particles. Dm < 100/um, 
decreased with increasing IWC values. For IWCs of order 0.001 g m“^, this portion from 
small particle was 80%), whereas for IWCs of order 0.1 g this fraction decreased to less 
than 30%. 

We fitted the PSDs generated by the MH97 parameterization to gamma distributions 
using Eqs. (4) to (7) to permit a qualitative comparison with our observations. The range 
of temperatures and IWCs at which the CEPEX and TRMM observations overlapped were 
examined: — 20 to — 50 C and 0.01 to 1 g m“^. Throughout this range single gamma dis- 
tributions did a good job of representing the two curves used for the MH97 parameteriza- 
tion, reliably capturing the inflections in the curves. 

It was difficult to compare the two sets of observations directly, because in MH97 the 
parameterization was derived in terms of the melted equivalent diameter which involved 
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assumptions about the particle density that differed from those made in this paper. Nev- 
ertheless, drawing upon the results of the comparison of A and the melted equivalent A in 
Fig. 19, we were able to note similarities and differences in p and Ap . Whereas our ob- 
servations indicated that Ap and p were primarily functions of T, the same coefficients 
using the MH97 parameterization were primarily functions of IWC and only somewhat 
functions of T. The range of /i were comparable, -1 to 2.5, for the fits to the MH97 param- 
eterization, whereas the Ap values from the parameterization were larger than would be 
expected by “melting” our spectra. More direct comparisons of the CEPEX and TRMM 
results will be the subject of a future study. 

6. SUMMARY AND CONCLUSIONS 

In this study, we have examined the characteristics and evolution of the PSDs in strati- 
form precipitating regions and anvils dmdng the TRMM field campaigns in Florida, Brazil 
and Kwajalein, Marshall Islands, emphasizing data from Kwajalein because the micro- 
physical probes used for those observations performed better. 

Cloud layers were sampled through the use of slow Lagrangian spiral descents, most of 
which commenced near cloud top and ended near cloud base (or at temperatures above 
0°C). This flight pattern facilitated investigation of how the PSDs changed in the vertical, 
both within the ice regions and in the melting layer. The patterns provided information 
on the vertical distribution of cloud microphysical properties, although the profiles could 
have changed temporally in a way partially characterized using radar data over the course 
of the spirals. While our Lagrangian spiral descents are far from perfect in assessing the 
evolution of particles in the vertical and in characterizing the properties of particle size 
distributions because of horizontal variability and temporal evolution, we believe that the 
results are a step towards increasing the understanding of how tropical ice cloud proper- 
ties vary in the vertical. 

The use of new instrumentation, allowing more accurate measurements of the concen- 
trations of larger particles and better definition of particle habits, has provided a data set 
that is more complete than has been previously available. The following summarizes the 
main findings. 

1. Relatively large particles were measured at cloud top, from 2 mm (990911) to 6 mm (990819, 
990830), even at T values as low as — 50° C. Such large particles have not been observed 
near the tops of layer clouds in the past. Particles of these sizes fall rapidly, above 1 m 
s“^. Rapid fallout, coupled with the particles’ increased growth downward through ag- 
gregation, which will generally further increase their fallspeeds, would indicate that there 
is rapid transport of moisture from upper tropospheric levels to mid and lower levels. 
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2. With distance below cloud top, the PSDs broadened, with some particles in the upper 
part of the ML reaching sizes exceeding 1 cm. Aggregation, which produced this growth, 
also led to depletion of particles smaller 0.1 cm, conforming to the observations of Field 
(2000) and others of how the aggregation process works. 

3. Changes in the characteristics of the PSDs were quantified through the use of gamma 

distributions of the form N — fitted to the PSDs, using a technique that 

matches the errors in the first, second, and sixth moments of the distributions. Expo- 
nential (p =0) were also htted to the PSDs. From cloud top to cloud base, or the top 
of the ML, the iVor and Ar values were fou n d to decrease monotonically. A similar ten- 
dency was observed for the exponential fits, with the results showing remarkable similar- 
ity to results from earlier studies of deep ice clouds as well as of cirrus clouds. Horizontal 
inhomogeneities over the course of individual loops of the spirals is manifested primar- 
ily in variations in A^or or A^o- The changed systematically with height, from positive 
(subexponential) in the upper parts of the clouds to negative (superexponential) in the 
lower parts of the clouds. The A^or and p values were highly correlated, and an expres- 
sion developed between the two parameters eliminated one unknown in the gamma dis- 
tribution form. 

4. The A^or values first decreased by about one order of magnitude near the top of the ML, 
whereas the Ap values decreased only slightly. Both the Ap and Aor values then increased 
abruptly in the lower part of the ML until the base was reached. 

5. The Ap and A values were generally a function of T, and curves were fitted to the data, 
although there was considerable scatter. The curve was similar to that reported by Ryan 
(2000) for mid-latitude synoptically and frontally generated ice clouds. We showed that 
the Ap versus T relationship could be used to estimate Aqp when either IWC or dBZg 
can be measured, although errors in Aqp of up to a factor of 3 might result from this 
method. However, this extent of error is no worse than errors resulting from the param- 
eterizations of the PSDs used in other studies. More accurate estimates may be possible 
through the use of vertically pointing Doppler radar data. 

6. Exponential curves fitted to our PSDs for tropical clouds that formed in association with 
convective forcing had A and Ng values that were similar to earlier studies for mid-latitude 
clouds formed primarily through synoptic forcing. The monotonic relationship found be- 
tween Nq and A or Agp and A may be applicable to a wide variety of clond types and ge- 
ographical locations. 

7. We developed general expressions between a number of macroscale and remote-sensing 
variables (e.g., IWC, Z, Vm, Vz, R, e, re) and Aqp and Ap . These expressions reduce 
to exponential if p is taken to be zero. The free variable in these expressions is Ag , a 
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term that will largely compensate for variations in IWC or dBZf, at a given tempera- 
ture. Terms needed to derive the coefficients in the relationship — specifically the depen- 
dence of ice particle density on size — were derived nsing some new techniques. Because 
these expressions are general, they can be modified as new information on ice particle 
density becomes available. 

8. It was shown that the variable Vz is fnndamentally a function of the spectral slope A 
or Ar- The distribution of Yz with dBZg we examined in the stratiform regions in Ap- 
pendix B provided implicit information on the distribution of A and Ap with dBZg. A 
similar analysis of ARMAR radar observations from three days from the TOGA COARE 
cirrus cloud observations, two of which appear in Heymsfield et al. (1998), show a simi- 
lar Vz vs dBZg pattern, suggesting that the distribution of A with dBZg is similar be- 
tween the two data sets. 

9. Although a plethora of ice particle shapes were observed — extending all the way from 
complex rimed crystals and rimed aggregates to pristine particles — much of the vari- 
ability could be explained by examining the Doppler radar data and quantified through 
use of the area ratio parameter. Rimed ice particles were sampled near deep convective 
clouds with sustained convection, as evidenced by the presence of deep and horizontally 
extensive anvils. Small (sub- mm), pristine particles were observed in clouds formed by 
more transient convection, with cloud-top temperatures controlling the habits: at low 
values of T, cirrus-type crystals, in some instances evolved into the capped columns that 
usually reflect growth at warmer temperatimes; at high values for T, planar-type crys- 
tals and capped columns dominated. Most particles of the larger sizes were aggregates. 
Relationships between the ice particle area ratio and diameter were used to describe the 
changes in the ice particle shapes with height, temperature, and Ar . 

10. Our PSDs, after “melting” them to obtain equivalent water PSDs and then fitting expo- 
nential curves to them, compared reasonably well with those of Sekhon and Srivastava 
(1970). We surmised that any differences resulted from the few samples (only 14) that 
were available to those researchers to captime the A versus dependence in the criti- 
cal regime where A values asymptote to about 9 cm“^. Also, about half of their samples 
were collected in the ML, where no direct relationship is evident between A and Nq val- 
ues. Although we could not directly compare our gamma size distributions to the tropi- 
cal cirrus PSD parameterization represented in terms of the melted equivalent diameter 
for T < — 20°C by McFarquhar and Heymsfield (1997), a cursory examination showed 
that these have similar characteristics. The latter parameterization is clearly superior in 
low-IWC regimes and at low temperatures, whereas the parameterizations we developed 
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are more flexible, are based on better measurements of the size distributions in larger 
sizes, and are better suited to the warmer temperatures and higher- IWC regimes. 

11 Given the range of conditions sampled, the results of our study apply to physically thick 
stratiform precipitation regions and anvils with temperatures from 3°C to about — 40°C. 
The observations apply to clouds falling into the category of “Deep Convection” accord- 
ing to the ISCCP optical-depth and height related cloud classification scheme. The ob- 
servations also apply to the range of IWCs from 0.01 to 1 g m“^, and radar reflectivities, 
as measured from overflying aircraft, in the range of 0 to 25 or 30 dBZe- 
The expression that we developed between IWC and the properties of the PSD (Eq. [9]) 
could have been substantially improved if direct measurement of IWC had been available. 
With direct measurements, assumptions that had to be made about ice particle density 
(the coefficients k, n, and a, and their variation within the ML) might not have been nec- 
essary. Also, questions about whetherAr depends on both temperature and IWC could 
have been addressed. Future direct measurements of IWC — along with PSD measure- 
ments from the probes — are needed. 

In closing, we would like to point out that there was considerable horizontal variability 
throughout most of our Lagrangian spiral descents and that the melting layers associated 
with four of our spirals were not horizontally extensive or homogeneous as might be the 
case in tropical mesoscale convective systems. A comparison of the ARMAR Vz versus 
dBZe relationships from the KWAJEX and TOGA COARE data sets suggests that the 
observations apply to a broad range of tropical ice clouds. However, more tropical obser- 
vations of PSD properties are clearly needed. We also want to point out that the TRMM 
satellite radar, which detects a minimum reflectivity of about 18 dBZg, would have seen 
some of the ice clouds we sampled on 990905 (TEFLUN-B), 990217 (LBA), 990819, and 
990830 (KWAJEX). Our observations may therefore be useful in understanding the prop- 
erties of the ice clouds and associated melting layers that the TRMM radar often detects. 
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Appendix A: List of Symbols 


Symbol Description 


A 

Ac 

A<P 

a, b 
af,bf 


C 

Co, Cl 
D 


Dm 

Dmm 

Dmax 

F 

9 

iwc 

k 

M 

m 

N 

Ni 

^0, Nor 
Nt 
n 
P 

PSD 


P 

R 

re 


Particle cross-sectional area; e.g., the particle’s area projected normal to its 
fall 

Total cross-sectional area of particle population (cm ) 

Area ratio— particle area divided by area of circle with the same D 
Coefficient, exponent in the fit of the area ratio versus diameter data from 
1-km imaging probe data 

Coefficients, exponents in Mitchell (1996) and H02 
Best Number versus Reynolds Number relationships 
Coefficient in the terminal velocity relationship 
Coefficient, exponent in equation relating X to T 

Maximum particle dimension, found by projecting its cross-section onto a 

horizontal plane normal to its fall direction 
Melted equivalent diameter 

Median mass diameter 

Maximum measured particle size in a given 10-sec sample 
Parameter related to the gamma fitting function 
Gravitational acceleration 
Ice water content 

Coefficient in equation for effective density 
Moment number for fitting function 
Ice particle mass 

Concentration per unit diameter as a function of D (cm 
Concentration in probe size bin I 

Concentration intercept parameter, exponential, gamma distribution 

_3 

Total concentration (cm ) of size distribution 
Exponent in effective density relationship 
Pressure (hPa) 

Particle size distribution 

Parameter giving the moment number for gamma fit 
Precipitation rate (mm hr ^ 

Radius used to describe the radiative properties of a particle population 
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r 

T 

VuVm.Vz 

y 

Z, dBZg 

a 

€ 

y 

K 

X, Ar 
y 

V 

Pa 

Pe 

Pi 


Correlation coefficient 
Temperature (°C) 

Ice particle terminal velocity, mean mass-weighted, and mean reflectivity- 

weighted velocities of ice paxticle size spectra 
Coefficient in mass versus diameter relationship 

Radar reflectivity factor, equivalent radar reflectivity (db) 

Exponent in effective density relationship 

Extinction coefficient 

Exponent in mass versus diameter relationship 

Exponent in the terminal velocity versus diameter relationship 

Slope parameter of size distribution for exponential, gamma PSDs 

Dispersion of gamma particle size distribution 

Kinematic viscosity 

Air density 

Effective ice particle density 
Density of solid ice 
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Appendix B: Estimation of Ice Particle Mass and Terminal Velocity 


The ice particle mass — a necessary input into calculations of the IWC, dBZe, and the 
mass-and reflectivity- weighted fall velocities {Vm and Vz) — can be derived from Eq. (1). 

As no measurements of the IWC were obtained to provide constraints on the estimated 
in Eq. (2), a comprehensive effort was undertaken to develop reliable estimates of pe and 
hence m. This effort involved three independent approaches: 

1. We compared calculations of the ter m i n al velocities with observations of snow crystal Vt 
values at the ground, a meaningful test of the pe relationship because Vt esc p°-^ to Pg-^. 

2. We compared the comprehensive set of dBZg and Vz values from the ARMAR (nadir 
viewing) on the DC-8 with calculations of these variables from the PSDs on the same 
days. The DC-8 generally avoided convection and measurements in convective cells that 
might have confounded this comparison. 

3. We compared calculations of rain PSDs, nsing as input snow PSDs near the top of the 
melting layer, with measured rain PSDs at the base of the ML during spirals. This method 
proved inconclusive, because of continued particle growth in the larger sizes and deple- 
tion in the smaller sizes. 

Our goal was to find a single relationship between pe and D across the full range of ob- 
served particle sizes to simplify the development of parameterizations of the PSDs in terms 
of macroscale variables (e.g., IWC) as presented in section 5. Three techniques (shown in 
Table 4) satisfied this requirement. The one parameter approach (pe cx; D) (Table 4), em- 
ploying the coefficients presented in Brown and Francis (1994) (hereafter BF94) which fit 
direct measurements of IWC and applied to ice clouds when particles were relatively dense 
and up to several millimeters in size.. We also tested the two-parameter “area ratio” (A^) 
approach, in which pe is a function of D and Ar (H02) (hereafter H02). It is our belief — 
and supported by the observations by Hanesch (1999) — that there is inherent information 
on Pe in the Aj , . The A^ two-parameter technique often reduces to a one parameter tech- 
nique because Aj. is usually approximately related to D by a power law. The two sets of 
coefficients shown for the Aj. technique in Table 4 were developed for low pe (Ay.[I]), and 
for relatively high pe associated with bullet rosettes up to about 0. 1-0.2 cm (Ay.[BR]). A 
third technique, a hybrid of the one-and two parameter techniques, was also attempted, 
as it captured the inherent information available in the two-parameter technique and the 
general decrease in pe with D, which is commonly found in transitions from single crystals 
at small sizes to aggregates at large sizes. For the hybrid technique, we used two sets of 
coefficients, consistent with earlier observations of aggregates at the ground (designated as 
H[I] and H[II] in the table. 
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From the images of some of the particles shown in Figs. 8-10 as well as and others not 
shown, some inferences about the may be made. Drawing upon observations in HOI, 
the hexagonal plate and platelike crystals in the 400-600-/rm range shown in Fig. 9. and 
larger than 800 pm in Fig. 15, have Pe values of 0.03-0.05 g cm“^. The capped columns 
in Figs. 8-10 have pe ~ 0.15 g cm“^, and the columns in Fig. 10 have pe ~ 0.07 g cm“^. 
Inspection of the particles smaller than 100 pm indicate that they are not solid ice spheres. 
The Pe for the more complex, partially rimed crystal shapes cannot be ascertained from 
the images alone. 

The Pe values as a function of D for the three techniques are plotted in Fig. 20A, using 
the coefficients given in Table 4. Using the second and third approaches, we derived Aj. 
from Ar — 0.4D“‘^ °®, which represents an average relationship for the TRMM spirals. 

Use of the BF94 approach produced unrealistically large pe values below about 500 pm 
and too small Pe values above 0.5 cm. Use of the ^1^(1) technique resulted in intermedi- 
ate Pe values at small sizes and relatively large pe values at large sizes. The Ar(BR) tech- 
nique, plotted over the applicable range of sizes, yielded relatively large pg. Using the two 
sets of coefficients (I) and (II) for the hybrid approach shown in Table 4 produced compa- 
rable results, with relatively large pe at small sizes and low pe at large sizes. 

The Vt was calculated as a function of D using the Pe from each approach. These calcu- 
lated terminal velocities were then compared to a large set of Vt measurements for snowflakes 
that fell to the surface in Japan at T of — 5°C (reported in Muramoto et ai, 1993) (Fig. 
20B). The shaded region in the figure shows the ± Icr bounds of the Vt observed in that 
study; are consistent with observations by many other We particularly plotted the data 
from the Muramoto study because these extend to somewhat smaller sizes. The Vt were 
calculated from Eq. (11). 

The calculated Vt generally increased with D (Fig. 20B), with less of an increase above 
about 0.4 cm, marking a change in the af and bf (drag-related) coefficients. The curve 
shown in the hgure for BF95 data suggests that values in that study were pe are overes- 
timated for aggregates in the small sizes and underestimated in the large sizes. The Vt 
trends with D for the A^ technique suggest that the pe values were are overestimated for 
D smaller than 0.05 cm using Ar(BR) and overestimated for D>0.5 cm for AR(I). The Vt 
trends with D for both hybrid approaches agreed well with the measured Vt for the larger 
sizes. The Vt below about 0.1 cm for the H(II) technique appear to be larger than observed 
for aggregates. 

As a means of evaluating the estimates of m from Eqs. (1) and (2), we used the PSDs 
for all spirals and Eq. (1) to calculate the dBZe (oc nrV) and Vz (related both to dBZg 
and terminal velocity) , then compared the results to the same parameters measured by 
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the ARMAR (nadir viewing) on the same days (Fig. 21). In the calculations of dBZe, 7.2 
dB was subtracted from the reflectivities to account for the differences in the dielectric 
constants between water and ice. Note that there may be some residual mesoscale ascent 
of order 10 cm s“^ just above the ML in these plots that is not averaged out. The mean 
values from ARMAR were obtained above the ML during those periods when the DC-8 
was in level flight and in most instances avoided convection. The four panels in the figure 
show the values calculated using the three techniques for obtaining p^. the single parame- 
ter technique pe oc here 6 is a power (n = 0 in Eq. [2]); the Aj. technique, pe = ^-4” 
[a — 0 in Eq. [2]); and the hybrid technique with two sets of coefficients. (The explicit 
equations used to derive Pe from the average TRMM versus D relationship are given 
in Table 4.) Although the standard deviations of the Vz are fairly large, a systematic in- 
crease of Vz with dBZ is noted, but with some natural variability and some variability re- 
sulting from vertical motions. 

Comparison of calculated and measured dBZ versus Vz indicates that when the parti- 
cles were small and the dBZ and Vz were relatively low, the one-parameter technique, us- 
ing the coefficients from BF95, produces results that overestimated Vz\ the Ar technique 
produced results that agreed moderately well with the radar data; and the hybrid (I) tech- 
nique’s calculations agreed well with the data for both sets of coefficients. When relatively 
large particles were present, at large dBZ and Vz, the one-parameter technique signifi- 
cantly underpredicted the Vz', the A,. (I) technique significantly overpredicted the Vz; and 
the two hybrid approaches produced reasonable Vz (I) and overpredicted the Vz (II). 

From these comparisons, we conclude that the hybrid (I) technique appears to give the 
most reasonable results. 
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Table 1 


INSTRUMENTS USED IN EARLIER TROPICAL ICE CLOUD STUDIES 


Study 

Probe 

Range {flm) 

Probe 

Range (flni) 

Probe 

Range (/im) 

Griffith et al. (1980) 

ID-C PMS^ probe 

70—460 





Knollenberg et al. (1982) 

2D grey 

40—2500 





Knollenberg et al. (1993) 

FSSP^ 

0.1-78 

2D grey 

50—1500 



Takahashi and Kuhara (1993) 

videosonde 

200- > 2000 





Pueschel et al. (1995) 

2D grey 

25-1800 





Heymsfield and McFarquhar (1996) 

FSSP 

0.1-21 

VIPS^ 

10-200 

2D— C 

40-1000 

McFarquhar and Heymsfield (1996) 
McFarquhar and Heymsfield (1997) 
Stith et al. (2001) 

FSSP 

2-45 ^m 

2D— C 

40-1000 

cpd 

<2 mm 


^ Particle Measuring Systems 
^ forward scattering spectrometer probe 

3 

video ice particle sampler 

4 

cloud particle imager 




Table 2 


TRMM Lagrangian Spiral Descents 


Date 

Location 

Times 

(UTC) 

Alt 

(m) 

Temp 
(° C) 

d{Alt) 

dt 

(m 

Loops 

(#) 

Op. Depth 
(Average) 

Radar 

Sfc. Precip 
(Any-]>18 dBZ) 

Change 

(dB) 

980905 

Fla^ 

213251 215040 

7923 3632 

-19 

6 

2.7 

6 

N/A 

Yes 

Yes 

< dz5 

990217 

Brazil 

194224 195601 

10042 4602 

-38 

0 

6.7 

4 

N/A 

Yes 

Yes 

< dz5 

990219 

Brazil^ 

201335 201934 

6091 4379 

-10 

2 

4.8 

3 

N/A 

Yes 

Yes 

ftilO 

990818 

Kwajalein 

035139 040701 

8508 5695 

-26 

-5 

3.1 

7 

24 

Yes 

No 

+5 

990819 

Kwajalein 

220448 225000 

6950 3370 

-14 

7 

1.3 

23 

37 

Yes 

Early 

— 25 

990822 

Kwajalein 

212049 220451 

11206 6983 

-50 

-16 

1.6 

14 

22 

No 

No 

Not detect. 

990823 

Kwajalein, 

031450 034550 

10341 6094 

-42 

-9 

2.3 

10 

25 

Yes 

Yes 

< ±5 

990830 

Kwajalein 

201056 203730 

7376 3685 

-18 

6 

2.3 

13 

29 

Yes 

Yes 

— 10 



205332 212120 

7319 3685 

-17 

5 

2.2 

13 

25 

Yes 

Yes 

— 25 

990911 

Kwajalein 

194955 203121 

10055 4514 

-39 

7 

2.2 

11 

30 

Yes 

Yes 

— 20 


1 2 

Near Cape Kennedy Near Puerto Velho N/A-Not available or unknown 



Table 3 

Equations to Derive Nor, Ar from Microphysical and Radar Variables 


Variable 

Full Equation 

Simplified Equation 

iwc [^] 

7Vor7rfca”xl0^r(4+fen+a+/A) 

3.3 X lo3jVorr(3.2+/i) 

QXr 

_^^(3.2+^) 

^mm 

3 + a + 6n+// + 0.67 

2.90 + a 

Ar 

Ar 

7 rnm£i 

Nav a^"xlO^^ r{7+2(6n+a) +/x) 

3.90JVor X 10’'r(5.5+/i) 

^ 1 m3 i 

p^2;)^j,(7+2(b«+c.)+/i) 

X^(5.5 + p) 

dBZe^ 

10 logioZ - 7.2 

10 logioZ - 7.2 

V™ [fi] 

CAr^”'^^r{K4'4+i*n+a!+/x) 

r(4+6n-f a+M) 

274r(3.8+(i)Ar~'’-^^ 

r(3.8+/i) 

[fl] 

CAr'~''>r(K+7+2(fcn+a)+/!i) 

218r(5.6+/i)Ar"‘’”‘^ 

r(74'2(6n+a)+/i) 

r(5.5+M) 


0.036/Vor7rfca” Cxl0®r(/c+44-6n+a:4-/i) 

3.2 X 10“iVorr(3.8+M) 

g]^^7(«+4+67r+aT7»l 

^p(3.8+ti) 

[^] 

iVo r 7f 1 0 ® r ( 3 + 6+ M ) 

4;^p(3+*>+M) 

2.3xlO®A^nrr(2.8+M) 

;Vj,(2.8+f«) 

e [km“^] 

0.2Ae 

0.2Ae 

re [/um] 

2^/{3)fert”“^xl0'*r(44-6n4-a4-/i) 

9ir(3.2-|-M) 


r(2.8+M)Ar«n 


^ Assumes Rayleigh scatters, valid for radar wavelengths of 2.1 cm and above. Units unless otherwise 
stated: cgs. 

Equations based on the following relationships: 

Mass: m = 

Effective density: = k{Ay)^D°‘, k = 0.04, n = 1.5, a = -0.5 

Area Ratio TRMM average values: a = 0.29, b = -0.18. 

Terminal Velocity: Vt = - 26/)2j[(3-l-6(n-i)+a)6/-i]^ 

Drag Coefficient Parameters: For Vm and R, o/ = 0.2072, bf = 0.638, for Vz, o/ = 1.0865, bf = 0.499. 
Kinematic viscosity: i/, Air density: pa 
Simplified equations for P = 500 hPa, T = 0°C 

K = (3 + b(n — 1) + a)bf — 1 

Fitted Functions 

Fitted Parameters Fit to Median Fit to median -la values Fit to median -|-la values 


Ar vs T(°C) 
<-27C 
A vs T 
p vs No 
p vs Ar 

Area Ratio=aD^ 


Maximum Diam. 


Ar = 6.8exp(-o-‘>®®^) 
Ar = 24.0exp(-o-O49D 
A = 12.5exp(-‘’o®®^) 
p = 0.21 ln(0.127Vor) 
p = 0.13Ar^° ®‘‘^ - 2 

a =0.32Ar^°°'‘* 
b =-0.45-F0.90o 

U„,,,=9.5A-o-9i 

D„,„,=3.2Ar-°-®^ 


4.0exp(-°-i°^) 

10.8exp(-”‘>®3T) 

9.3exp(-°°®3U 

p = 0.21 ln(0.04Aor) 
p = 0.14Ar^“ ®®' - 2 

0.31Ar 

-0.42-l-0.75a 

8.3A-'>»4 

3.4Ar 


11.5exp(^°'°®®^) 
53.3exp(^°'°^^^) 
16.7exp(^°'°®®^) 
p = 0.21 ln(0.33Vor) 

p = 0.11Ar^°'^^^ - 2 

0.33Ar^°°^® 

-0A7+lMa 



Incomplete gamma function: •y{a,x) = /q^ 

Example: IWC = + bn +a.+M,ArC„„.) 







Table 4 


Parameters Used to Calculate pg 


Technique 

Equation 

Designation 

Coefficients 


Simplified Eqn. 

One-Parameter 

pe = aD^ 


a = 

^ .00561 

II 



Two-Parameter 

pg = k{Ar)'^ 









Ar(I) 

k 

^ 0.18 

n ^ 2.50 

Pe 




Ar(BR) 

k 

^ 0.48 

n ^ 2.25 

Pe 


Hybrid 

Pg = k{ArTD<^ 









H(I) 

k 

^ 0.03 

n ^ 1.2 

Pe 

= 0.0067U^°-®°* 




a 

^ -0.4 






H(II) 

k 

- 0.04 

II 

e 

Pe 





a 

^ -0.5 





See Appendix B for discussion of nomenclature 



Figure Captions 


Figure 1. Flight tracks of DC-8 and Citation during Lagrangian spiral descent on 0819 (top) and 0823 ( bottom ), 
both from KWAJEX. Period of DC-8 overpass of Citation track is shown with bold line. 

Figure 2. Radar data as a function of height from ARMAR (top five rows) and EDOP (bottom two rows) av- 
eraged over time periods when DC-8 or ER-2 crossed over or in proximity to Citation track. Radars are nadir 
viewing and times are indicated in left panels. Left; Mean measured dBZg, with zbl (T bounds over periods 
shown. Right panels: Doppler velocity, and zhl(T bounds of measured Vz over periods are shown. 

Figure 3. Average PSDs for various loops of the 990819 spiral. Heavy lines : 2D-C. Light lines : HVPS. 

Figure 4. Measured reflectivities ( left ) and representations of the PSDs ( right ) during Lagrangian spiral de- 
scents on three days during KWAJEX. ARMAR measurements were obtained from DC-8 in nadir-viewing 
mode, with abscissa showing time. The boxes ( left panels ) show the top, bottom, and approximate horizon- 
tal positions of the Citation as if the spiral occurred during the overpass period (see text). 

Figure 5. Same as Fig. 4, except that bottom panels shows data from EDOP on the ER-2 (data courtesy of 
Gerry Heymsfield) during LBA. Data for 990830 is from the first spiral listed in Table 2. 

Figure 6. Examples of particles imaged by the HVPS probe for the 14 loops (numbered in panels) of the spi- 
ral on 990822. The scale is shown in the top left panel. 

Figure 7. Concentration of particles 100 fJ, and above as measured by the 2D-C probe during five of the La- 
grangian spiral descents. The distinct spikes at some points in the TVj’ traces in the ML are probably an arti- 
fact of the 2D-C probe. 

Figure 8. Examples of particles imaged in three size ranges by CPI probe on 990819. Magnification between 
different size ranges is not the same. 

Figure 9. Same as Fig. 8, except for 990822 (KWAJEX). 

Figure 10. Same as Fig. 8, except for 990217 (LBA, top ) and 0823 (KWAJEX, bottom ). The 990217 case 
shows riming near to convection and the 990823 case shows pristine cirrus-type crystals reflecting short-lived 
convective elements and low temperatures 

Figure 11. Microphysical variables derived from PSDs for three spirals. Rows from top to bottom : IWC, R, 
dBZg as calculated assuming equivalent water spheres, Vzi ^■nd V^. 

Figure 12. Fit parameters for exponential ( left ) and gamma ( right ) distributions as functions of T from 1- 
km PSDs. Top panels : A and Ap. Fit from Ryan (1996) and to the TRMM data in left panel, and to TRMM 
data in right panel, shown with solid line. Fits exclude 0822 case. Middle panels : Nq and Nqy- Bottom left : 



correlation coefBcient. Bottom righ t: /i versus T . Different colors represent different days. Inserts in panels 
B, C, D: data for 0911 case. 


Figure 13. Findings from the gamma distribution fits. A: Exponential versus gamma A values. Solid line, 1:1 
line; dashed: Fit. B: Ap calculated from T versus measured. C, D : jJ, versus Ap, fJ- versus iVpp. Mean and 
ih 1(7 bounds in equal intervals of abscissa variable. Fits to the data (excluding 0822 case). 

Figure 14. Maximum particle diameter per 1 km (about 8-9 sec) PSD as a function of A ( top i and Ap ( bottom ). 

Figure 15. Nq versus A ( top ) and iVpp versus Ap for the various spirals. Large symbols represent data points 
at 0 and 3°G. 

Figure 16. Coefficients in the relationship Aj- = 0,D^ as derived by fitting the Aj- versus D data from 

the imaging probes. A, B: Height dependence. C, D: Temperature dependence. E: Dependence of O on A. F: 
Relationship between CL and b coefficients. 

Figure 17. Ratios of various parameters as obtained from the fitted size distributions to those calculated from 
the PSDs. 

Figure 18. The Nq versus A points measured during Lagrangian spirals by Gordon and Marwitz (1986), Lo 
and Passarelli (1982), Passarelli (1978), and Field (1999). Unless otherwise noted, data are from descents. The 
listed values in GM86 (2D-C probe) and the values extracted from the graphical plots in P78 and LP82 (ID-P 
probe) are plotted. Temperatures at top and bottom of each spiral are shown. The dotted lines show IWCs 
calculated using the Brown and Francis (1995) pe~D relationship. The dashed lines in the lower right panel 
show an envelope of the observations from TRMM. The lines labeled “Ryan” use the A-T relationship from 
Eq. (8) to predict A from T, plotted at the measured value of Nq. 

Figure 19. Comparison of the Sekhon and Srivastava (1970) parameterizaton (solid line) with the TRMM re- 
sults (points). The line was obtained by the method described in section 5, with the squares along the line 
corresponding to the four precipitation rates that were used to develop the parameterization. The unmelted 
Nq — — A points and the ’’melted equivalent” points were obtained from the 990822 spiral descent. 

Figure 20. Top : The usingthe one-and two-parameter techniques as well as the hybrid technique. Bottom : 
Corresponding V^, compared with data from Muramoto et al. (1993), shaded region. 

Figure 21. Comparison of measured V^-dBZg from ARMAR radar (mean and ihl (T bounds) with values 
calculated from PSDs (assuming equivalent water spheres) using (A) one-parameter technique with Brown and 
Francis (1995) coefficients, (B) Aj. technique, (G) H(I) technique, and (D) H(II) technique (coefficients listed 
in Table 4). 
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A— Dependence, Lagrangian Spirals in Larlier Studies 


GM86, 15 Feb. 1982 


GM86, 13 March 1984 



10 

A (cm”') 


10 

A (cm”') 


LP82, 25 Feb. 1980 


LP82, 26 Feb. 1980 


F-6' ^,-29, 


N O/ 

.O'. . .O'.-, . . 

10 

A (cm”') 

LP82, 8 March 1980 


N or 
.O'. .-. . .O'.-, <y.'. 

10 

A (cm”') 


P78, 10 March 1975 (ascent) 


-A p; 
or 

■X G-’ 

O' ^or 


10 

A (cm”') 

P78, 10 March 1975 


§ 0.100 



or ^r 
.O'. .-. . .O'.r -y.'. 

10 

A (cm”') 


F99, 25 October 1996 (spiral C) 



TRMM / 
'f 

/ r' 


/ / /) 7.' 

/ r' / J^9 G' 




mil 










Comparison with Sekhon and Srlvastava Parameterization 













